The Cosmology of Modified Gauss-Bonnet Gravity 
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We consider the cosmology where some function f{G) of the Gauss-Bonnet term G is added to the 
gravitational action to account for the late-time accelerating expansion of the universe. The covariant 
and gauge invariant perturbation equations are derived with a method which could also be applied to 
general f{R, R^^Rab, R"'"^'^ Rabd) gravitational theories. It is pointed out that, despite their fourth- 
order character, such f{G) gravity models generally cannot reproduce arbitrary background cosmic 
evolutions; for example, the standard ACDM paradigm with Ode = 0.76 cannot be realized in f{G) 
gravity theories unless / is a true cosmological constant because it imposes exclusionary constraints 
on the form of /(G). We analyze the perturbation equations and find that, as in f{R) model, the 
stability of early-time perturbation growth puts some constraints on the functional form of f{G), in 
this case f /dG^ < 0. Furthermore, the stability of small-scale perturbations also requires that / 
not deviate significantly from a constant. These analyses are illustrated by numerically propagating 
the perturbation equations with a specific model reproducing a representative ACDM cosmic history. 
Our results show how the /(G) models are highly constrained by cosmological data. 
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I. INTRODUCTION 

The accumulating astronomical evidence for the 
present accelerating expansion of the universe has stim- 
ulated many investigations into the nature of the dark 
energy which might be responsible for this unexpected 
dynamics (for a review see, e.g., Jj). Besides proposing 
to add some kinds of exotic (and purely theoretical) mat- 
ter species into the energy budget of the universe, many 
investigators have also focused their attentions on modi- 
fying general relativity (GR) on the largest scales so as to 
introduce significant modifications at late times. One ex- 
ample is provided by the family of f{R) gravity models, 
which had also been considered before the discovery of 
cosmic acceleration (see for example Refs. [l, 0, Q)- In 
Refs. 0, the authors discuss a specific model where 
the correction to GR is a polynomial function of the 
B?,R'^^Rab and R'^'^'^'^Rabcd quadratic curvature invari- 
ants (here, R, Rab and Rated are respectively the Ricci 
scalar, Ricci tensor and Riemann tensor) and showed that 
there exist late-time accelerating attractors in Friedmann 
cosmological solutions to the theory. It is very interest- 
ing that when the Ricci scalar R in the Einstein-Hilbert 
action is replaced by some general functions of R and 
R'^''Rab, it becomes necessary to distinguish between two 
different variational approaches in deriving the field equa- 
tions. In the metric approach, as in Refs. ^, i6|], the met- 
ric components gat are the only variational variables and 
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the field equations are generally of fourth-order, which 
makes the theories phenomenologically richer but more 
stringently constrained in many cases. Within the Pala- 
tini variational approach, on the other hand, we treat 
the metric gat and connection F^^ as independent vari- 
ables and extremize the action with respect to both of 
them, and the resulting field equations are second order 
and easier to solve. The Palatini f{R) gravity is also 
proposed as an alternative to dark energy in a series of 
works 7, 8, 9, lOj. There has since been growing interest 
in these modified gravity theories: for the local tests of 
the Palatini and metric f{R) gravity models see [lllfl^. 
and for the cosmologies of these two classes of models see 

Both approaches to modifying gravity are far from 
problem-free. In the metric f{R) gravity models, the 
theory is conformally related to standard GR plus a self- 
interacting scalar field W\ , which generally introduces ex- 
tra forces inconsistent with solar system tests The 
Palatini approach, on the other hand, generally leads to 
a large (or even negative) sound-speed-squared term in 
the growth equation of the matter perturbations on small 
scales, and its predicted shapes of the cosmic microwave 
background (CMB) and matter power-spectra deviate 
unacceptably from those predicted in GR OIUjII^i and 
so fail the most fundamental cosmological tests. Again, 
these examples reiterate the difficulties encountered when 
trying to make general modified gravity theories compat- 
ible with observations. 

In this work we will focus on another form of mod- 
ified gravity, the modified Gauss-Bonnet theory, which 
is proposed and discussed in [10, HJ, |^ 0, 01 (see, 
for example, [H, [H, [23| for a related model). In these 
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models, GR is modified by adding an arbitrary func- 
tion /(G) into the gravitational action, where G = 
g2 4^Jiabji^^ _,_ R^bcdji^^^^ ^Yie Gauss-Bonnet invari- 
ant (which is a topological invariant in four dimensions). 
Such correction is motivated by the effective low-energy 
actions in string theory and is shown to be able to pass 
solar-system tests even though it arises in the metric vari- 
ational approach. Furthermore, refs. [13, [2^ also demon- 
strate that such models can produce late-time cosmic ac- 
celeration, as well as a transition from deceleration to 
acceleration, or from a non-phantom phase to a phan- 
tom phase. Here, we are interested in the perturbation 
dynamics in such modified gravity theories, and the cos- 
mology arising from them at first-order in perturbation 
theory. 

Our presentation is organized as follows. In Sec. [TTl 
we briefly introduce the model and present the pertur- 
bation equations of general /(G) model in covariant and 
gauge invariant (CGI) form. For this, we generalize the 
effective energy-momentum tensor approach of deriving 
CGI modified gravitational field equations in a way which 
could be apphed to a general f {R, R"'' Rat, R"'""^ Rated) 
theory - which includes the /(G) and /(i?) models as 
specific cases. In Sec. lIIII we shall discuss the background 
evolution of /(G) models and explain why arbitrary cos- 
mic histories cannot be realized with the /(G) model. 
The perturbation equations are then analyzed in Sec. lIVi 
where we also evolve these equations numerically to in- 
vestigate the /(G) effects on the growth of linear pertur- 
bations. Our discussion and conclusions are presented 
in Sec. [V] Throughout this work our convention is cho- 
sen as [Va, 'Vbju" = R^^ju'^, Rab = Racb where a,b,-- ■ 
run over 0, 1,2,3 and c = h = 1; the metric signature is 
(+,—,—,—) and the universe is assumed to be spatially 
flat and filled with photons, baryons, cold dark matter 
(CDM) and three species of effectively massless neutri- 



II. FIELD EQUATIONS IN MODIFIED 
GAUSS-BONNET GRAVITY 

In this section we briefly introduce the main ingredi- 
ents of /(G) gravity and derive the general perturba- 
tion equations that govern the dynamics of small inho- 
mogeneities in the cosmological models that arise in this 
theory. 



A. The Generalized Einstein Equations 

Our starting point for /(G) gravity is the modified 
Einstein-Hilbert action. 



S 



R + fjG) 
2k 



(1) 



action with respect to the metric gab gives the modified 
Einstein equations 

Rab — -^9abR 

= + \gabf - 2FRRab + AFRlRbc 

— 2FRacdeRb'^'^ — ^FRacdbR"^ + 2RV bF 

-2RgabV^F - 4i?^VfcV,F - 4i?gV„V,F 
+ARab^^F + AgabR'-'V^^dF - ARacbd^^V^F, (2) 

where F = F{G) = df{G)/dG (notice that, unhke in the 
f{R) models, F here is not dimensionless) and T^j, is the 
energy-momentum tensor of the fluid matter (photons, 
baryons, cold dark matter, and light neutrinos). The 
trace of Eq. ^ reads 



-R = K{pf -ipf) + 2f -2FG 
~2RDF + 4i?a6V° V''i^, 



(3) 



in which k = SttGn, with Gn being the gravitational 
constant and R = R{gab) is the Ricci scalar. Varying this 



with = T = — 3p^ . We see that the curvature- 
related quantities (R, G,F, ■ ■ ■) are determined by the 
energy-momentum tensor of the fluid matter through a 
complicated dynamical equation Eq. ^ and thus the 
modification to the GR field equations can be understood 
as a change in the way that the spacetime curvature, and 
thus the Ricci tensor Rab, responds to the distribution of 
matter. 



B. The Perturbation Equations in General 
Relativity 

The CGI perturbation equations in general theories 
of /(G) gravity are derived in this section using the 
method of 3 -I- 1 decomposition [H, [l^, [s^, [3l|. Fur- 
thermore, we shall adopt the effective energy-momentum 
tensor approach [l^, (3^ , which treats the modifications 
on the right hand side of Eq. ^ as an effective energy- 
momentum tensor. However, since the modification gen- 
erally involves terms nonlinear in Rab, Rabcd, we should 
express these terms appropriately. This will be done in 
more detail below, but now let us briefiy review the main 
ingredients of 3 -I- 1 decomposition and their application 
to standard general relativity for ease of later reference 
(as in GR, there is only fluid matter so we shall neglect 
the superscript ^ in this subsection). 

The main idea of 3-1- 1 decomposition is to make space- 
time splits of physical quantities with respect to the 4- 
velocity of an observer. The projection tensor hab is 
defined as hab = gab — UaUb and can be used to obtain 
covariant tensors perpendicular to u. For example, the 
covariant spatial derivative V of a tensor field Tj^'..^ is 
defined as 

sj-Tt:,i^Kh^--- hth:^--- KN'T}::1 (4) 

The energy-momentum tensor and covariant derivative 
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of the 4- velocity are decomposed respectively as 

Tab = T^ab + '2q{aUb) + pUaUb ~ phab, (5) 
WaUb = (7ab + T^ab + ^dhab + UaAb. (6) 

In the above, nab is the projected symraetric trace-free 
(PSTF) anisotropic stress, Qa the vector heat flux vec- 
tor, p the isotropic pressure, cTab the PSTF shear ten- 
sor, zuab = V[aUb], the vorticity, = V'^Uc = 3d/ a {a is 
the mean expansion scale factor) the expansion scalar, 
and Ab — iib the acceleration; the overdot denotes time 
derivative expressed as (j) — u°'V a4>, brackets mean an- 
tisymmetrisation, and parentheses symmetrization. The 
normalization is chosen to be it°Uo = 1. The quantities 
TTab, Qa, P,P are referred to as dynamical quantities and 
Cab, ra7ah, ^, as kinematical quantities. Note that the 
dynamical quantities can be obtained from the energy- 
momentum tensor Tab through the relations 



the Weyl tensor Wabcd, defined by Sab = u'^u'^^Vacbd and 
I3ab = —\W^u'^£ac^^£fbd- The angle bracket means tak- 
ing the trace-free part of a quantity. 

Besides the above equations, it is useful to express the 
projected Ricci scalar R into the hypersurfaces orthogo- 
nal to u° as 



^ 3 



(18) 



The spatial derivative of the projected Ricci scalar, rja 
^aVa-R, is then given as 



2a - 

and its propagation equation by 



(19) 



2,9 2 
??a + y??a = --eaVaV -A-anVa^-q. (20) 



p ^ TabU'^u", 

P = -Ih'^'Tab, 



TTab = h^hf Ted + phab- 



(7) 



Decomposing the Riemann tensor and making use the 
Einstein equations, we obtain, after linearization, five 
constraint equations [13, HH : 



= V%e''\^u''wab); 



KQa 
Bab 



^ h V CTafc + 

1 2 „ 



^b)ec'^^i 



2^ 



Vc^d + (p + p)n7cd 



cd^,b 



(8) 
(9) 
(10) 

(11) 
(12) 



and five propagation equations, 

O + l0^-^^Aa + ^{p + 3p) = 0;(13) 
2 - 1 

<Jab + -^Oaab - (aAb) + Sab + ^UTTab = 0;(14) 

2 

■dj + -em-WiaAb] = 0;(15) 



0;(16) 
0.(17) 



1 






1 




T^ab H 










"2" 



Sab + dSab — V'^Srf(ae^-)g 
Bab + OBab + ^"SMaS 



u 

d{at-b)ec " 
d 



d{a^b)ec 



U 



Here, Eabcd is the covariant permutation tensor. Sab and 
Bab are respectively the electric and magnetic parts of 



Finally, there are the conservation equations for the 
energy-momentum tensor: 

p+{p + p)e + V-'qa = 0, (21) 

<7a + ^fga + (P + p)Aa - VaP -I- VVah = 0. (22) 

As we are considering a spatially-flat universe, the spa- 
tial curvature must vanish on large scales and so i? = 0. 
Thus, from Eq. (fT8|l . we obtain 



\0^ = np. 



(23) 



This is the Friedmann equation in standard general rel- 
ativity, and the other background equations (the Ray- 
chaudhuri equation and the energy-conservation equa- 
tion) are obtained by taking the zero-order parts of 
Eqs. (HSIIIT]), as 



p + {p + p)e = 0. 



(24) 
(25) 



In what follows, we will only consider scalar modes of 
perturbations, for which the vorticity voab and magnetic 
part of Weyl tensor Bab are at most of second order [s^, 
31j and will be neglected from our first-order analysis. 



C. The Perturbation Equations in /(G) Gravity 

In the effective energy-momentum tensor approach, 
the field equations Eqs. ([S] - [25]) listed above preserve 
their forms, but the dynamical quantities p,p, qa, i^ab ap- 
pearing there should be replaced with the effective total 
ones p 



TT^jj -|- TT^j, in which a superscript means the contribu- 
tion from the Gauss-Bonnet correction. 
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Writing the modified Einstein equations, Eq. ([2]) in the 
following form, 



Rab 



abi 



(26) 



one can easily identify 

«7^a1 = ^gabf~2FRRab + 'iFRlRi,c 



-2FRacdeRb 



cde 



AFRacdbR"^ + 2RVaVbF 



-2Rgaby^F ~ ml^bVcF ~ AR^^a^cF 

+ARabV^F + 45afci?"'*VcVdF - ARacbdV^V^F 

and Eq. ([7]) can be used to calculate , p*^ , 5^ , tt^^ . In 
order to do this we need the explicit expressions for Rab 
and Rabcdi which could be obtained in terms of either (ef- 
fective total) dynamical quantities or kinematical quan- 
tities, or a mixture of the two. 

To express Rab and Rabcd explicitly, now decompose 
the symmetric Ricci tensor Rab in the following general 
way. 



Rab = AUaUb + =Lhab + 2u(a^b) 

then Eq. ((26)) gives 



(27) 



A = 



3p*°*) 



2 ^'^ 



'2 - V'^Aa 



1 



T — 



3 



ptot) 

+ R- V^Aa 



= -2 



f^ab 



(28) 



in which we have used Eqs. (HI [131 [H [HI) . Notice that 
the first lines are expressed in terms of total dynamical 
quantities and the second lines of kinematical quantities. 
For those terms involving Rabcd: we shall use the decom- 
position of Riemann tensor extensively (keeping in mind 
that u'^u'^Wacbd = £ab)- 

1 

Rabcd = -^[gacRbd + gbdRac — gadRbc — gbcRad) 

+yVabcd - -^Rigacgbd - gadgbc)- (29) 

6 

For example, it is easy to show that, up to first order, 

Ra'^'^Rbcde 



= —-^d^ab + -^dabR'^^Rcd + -^RRab 



7a6i?', (30) 



SO that 



^abcd — Kab — g-K , 



C I?2 o nab D 

Lt — —it — zrC tiab 



(31) 



where 



R = -26- -9'^ + 2V°-Aa~ R] 



U Kab 



V'^Aa. (32) 



With these useful relations and some calculation, the 
contribution to the energy-momentum tensor from the 
Gauss-Bonnet correction term can be identified as 



'^P'' - -if-FG) + -iA-3E)F6 



--iA-3E)V^F; 



-^P"" = -{f-FG) + -{A-3E)F 



A{9F + W^F); 
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-^(A-3S) \^aF-^e^aF 



(33) 



(34) 



-FOTaiih) 



-A 



[Paab 



-AFSab + -eFEab. 



(36) 



Here, we want to make some comments about these equa- 
tions. First, if / is constant, then F and its derivatives 
vanish, so that np^ = —Kp^ — ^ and Kq^ = ktt^^ = 0, 
and thus we have the ACDM limit. Second, it is not 
difficult to check that the above quantities satisfy the 
independent energy-momentum conservation equations 



p° + (pG+p«)^? + V"9a'' = 0, 

(pG+pG)Aa-VapG + V''7rS, = 0. 



This is a result of the energy-momentum conservation in 
fluid matter and the contracted Bianchi identity. Thirdly, 
it would be convenient to use — ^{ql + q^) and 
Sab = K{n[f^ + TT^J (see Eq. ^) to rewrite Eqs. jMl 
l36l) so that Kq2 , ktt^^ are expressed respectively in terms 
of Kql, K7r{(, and other quantities. This is what we do in 
the numerical calculation. Fourthly, it is interesting to 
note that there is no Va^ term in KVa/c*^ up to first order 
because A - 3S = f^^ _^ i? and Va(/ - FG) ^ -GVaF; 



this is positive because otherwise Eq. p9|) will no longer 
be an algebraic equation for Vq0. For similar reasons 
Eq. ([20]) remains a first-order differential equation for 
rja for the present model. These simplifications also 
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occur in the f{R) gravity models but not in general 
f{R,R''''Rab,R°-'""^Rabcd) theories. In the later case, the 
method we use here to derive the CGI perturbation equa- 
tions still applies and the perturbation equations will be- 
come even higher order and more complicated (specifi- 
cally, some of the perturbation equations above will be- 
come propagation equations for q^'^ and tt^*-^ , where 
denotes general modified gravity theory 33,iM])- Finally, 
we can see that the quantity F here appears only to (at 
least directly) influence background evolutions, and with 
the background fixed it is its derivative which determines 
the perturbation evolutions. 



III. THE BACKGROUND EVOLUTION IN /(G) 
MODELS 



In this section we discuss the background evolution in 
general /(G) gravity models. Recall that in f{R) gravity 
theories the fourth-order nature of the Friedmann equa- 
tion allows enough freedom for the model to reproduce 
an arbitrary background cosmic evolution. Since the 
field equations are also fourth order in the /(G) models, 
one might think that they could also describe arbitrar- 
ily parameterized background histories. However, this is 
not the case, as we shall see below. This is because G 
must change its sign in the recent past in many fixed- 
background models. 

For the background evolution, we use the Friedmann 
equation Eq. (I23|) with the dark energy density given by 
Eq. 



1, 



1 



-9^F (37) 



where we have kept terms only up to zero order and used 
the expression G — ^9^6 + SjO\ which is obtained from 
Eqs. dlZlllHlEIl)- Following we define the following 
dimensionless quantities (here H = 6/3 is the Hubble 
rate and Hq is its present-day value) 



E 



IP 

^0 



(38) 



in terms of which G, and so Eq. (j37p can be written as 



G = 12H^E{E* + 2E), (39) 



with 



f G G**\ , G* 
^l^DEexp [-3(l + w)iV], 



(40) 



where a star denotes the derivative with respect to = 
log(a), JIde = kpy)e/3Hq is the dark energy fractional 



energy density today, and w = const, is the usual dark- 
energy equation of state (EOS) parameter. Note that by 
writing in this way we have chosen to parameterize the 
background expansion to be the same as the dynamical 
dark-energy model {w ^ —1) or the ACDM paradigm 
(^^ = -1)- 

In the following, we shall assume w = —1 for simplicity 
and the calculations for general w could be done similarly. 
In this case we have 

E = Ode + f^mexp(-3iV) + f7^exp(-4A^), (41) 

where fi™, il^a are respectively the fractional densities for 
non-relativistic and relativistic matter species. Deep into 
the radiation-dominated era we have E = ftr exp(— 4iV), 
and thus Eq. (liO]) reduces to 



y** + 9y* + 8y 
whose solution is 



de, 



(42) 



y{N) = Aexp(-A^) + Bcxp(-8A^) + 6r2DE. (43) 



In this work we shall require 

/(G) 



lim 



R 







so that B could be set to zero in Eq. To obtain the 

background evolution numerically, we start deep in the 
radiation-dominated era {e.g., at a = 10^^) and take the 
radiation-dominated solution, Eq. (j43p . to be the initial 
condition. The solutions to Eq. are then character- 
ized by a single parameter A. Different values of A give 
the same background history, but in general lead to dif- 
ferent evolutions for the perturbations, as we shall see 
below. Note that A = describes the standard ACDM 
paradigm. 

Up to this point the procedure is quite similar to 
that in metric f{R) gravity models. However, as was 
claimed above, the /(G) gravity model cannot be used to 
reproduce arbitrarily parameterized background expan- 
sion histories. To explain why, we shall again take the 
ACDM background as an example and adopt the value 
^DE — 0.76, as suggested by the Wilkinson Microwave 
Anisotropy Probe (WMAP) three-year data (sH]. Then 
it will be easy to find that G* changes its sign (from -|- 
to — ) at A^o ~ —0.153, which means that G increases 
(decreases) when N < (>)7Vo- As a result, with the 
match to ACDM at N < Nq the function /(G) with 
all possible values of G has been determined, and there 
will generally be no freedom left to fix the evolution to 
ACDM at > iVo as well (of course, if /(G) is a real cos- 
mological constant then ACDM will also be reproduced 
at > A'o, but in general reproducing ACDM on both 
sides of Nq is far too strong a requirement to be satis- 
fied). So, what we may conclude is that the /(G) model 
can mimic a ACDM universe up to Nq, after which the 
evolution might be governed by the already-determined 
f{G). However, note that A'o ~ —0.153 corresponds to a 
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FIG. 1: The evolution versus redshift of f{G)/R for the /(G) 
gravity models fixed to match a ACDM cosmic expansion 
history {Qde ~ 0.66). The curves from top to bottom are 
characterized by ^ = -0.02,-0.01,0 (ACDM), 0.01, 0.02 re- 
spectively. Notice that with our convention R < and / > 0. 



critical redshift zq ~ 0.166, so the transition from ACDM 
phase to a non-ACDM phase occurs quite late [s^ . 

In the above, we have analyzed for a ACDM (w = — 1) 
background. Nonetheless we can expect that similar situ- 
ations exist for general dynamical (w ^ — 1 ) dark-energy 
backgrounds which are characterized by a late-time tran- 
sition from G* > to G* < 0. The investigation of 
background evolutions in general modified Gauss-Bonnet 
models is an interesting topic which is unfortunately be- 
yond the scope of the present work. Here, we solve the 
perturbed equations numerically and investigate the ef- 
fects of the /(G) modifications to GR on the growth 
of linear perturbations. For simplicity we shall adopt 
a slightly unrealistic ACDM cosmic history which is de- 
scribed by r^DE = 0.66 (in this case the transition from 
G* > to G* < has not taken place yet). 

In Fig.[Tl we show the redshift evolutions of /(G) com- 
pared with R. Each curve here is characterized by a 
specific value of the coefficient A, which can be either 
positive, negative or 0. Note that although all these 
curves lead to the same ACDM background evolution, 
not all of them are cosmologically viable. As will be 
discussed below, the stability of early-time perturbation 
growth requires Fq < (where Fq = dF/dG) and that 
of small-scale perturbation growth furthermore requires 
\Fg\H^ < 1. The condition Fg < is found to corre- 
spond to the A> Q subclass of the solutions to Eq. (j40| . 
The evolution of \Fg\H^ = ~FgH^ (which is dimension- 
less) in this subclass are shown in Fig. [2] It is obvious 
that \Fg\H'^ is a rapidly increasing quantity with respect 
to time whose magnitude is much smaller than 1 at early 
times. 



FIG. 2: The evolution versus redshift of the dimensionless 
quantity ~FgH^ for the /(G) gravity models fixed to match 
a ACDM cosmic expansion history (J^de = 0.66) which 
are stable in early time perturbation growths [Fg < 0). 
The curves from bottom to top are characterized hy A = 
0.0001,0.001,0.005,0.01,0.02 respectively. Note that a true 
cosmological constant has Fa ~ 0. 



IV. THE PERTURBATION EVOLUTION IN 
/(G) MODELS 

Now we turn to the evolution of perturbation equations 
in the /(G) model. The equations presented in Sec. Hi] 
are not yet closed and from Eqs. ([33] -[SG]) we see that 
an evolution equation for V^F is also needed. To obtain 
this, let us look at the trace equation, Eq. ^ (which can 
also be obtained from Eq. (19) and the spatial derivative 
of Eq. (13)). Taking its spatial covariant derivative, af- 
ter some manipulations we obtain the following evolution 
equation 



= e 




271 



iFe 



16 e^Fo 



(44) 



where e is the harmonic expansion coefficient of Vq-F, as 



(45) 



and Q'^ = f VqQ'^ and Q'' are the zero-order eigenval- 
ues of the comoving Laplacian a^V^ {a?V'^Q^ = k'^Q'^), 
Fg = dF/dG = 9^//(9G^, and the source function S is 
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given by 



S 



3k 
1^ 



{Xf - 3XPf) - 3 



1-Ifo 



6*3 a' 




(46) 



where , X^^ , Z,ri are respectively the harmonic ex- 
pansion coefficients of VaP^ and (see 
for example [3lj|). Moreover, here we are working in the 
CDM frame (with the 'observer' comoving with dark- 
matter particles and so free-falling) in which case we can 
set Aa = [13, [Ml to simplify computations. In this case 
we have (up to first order) 



^ — ^ n 



table. As we have seen in the above section, the magni- 
tude of the dimensionless quantity I^^-FgI is tiny at early 
times (deep into the matter- and radiation-dominated 
eras) so that at that time this term dominates over the 
other two in the squared brackets. If Fq > 0, then 

~TE''W^r ^ these early times, which makes 

the perturbation e unstable and it grows quickly to spoil 
the linear theory. This is similar to the analysis of sta- 
bility in f{R) models For the subclass of mod- 
els with Fq < 0, to which we are restricting ourselves, 

27 1— - F6 

~iE e-^fc ~* °^ value of e quickly settles to- 

wards — ^9'^FcS . It can be checked easily that this is 
equivalent to WaF = Fa'^aG, which is as expected (note 
that at these times the influence of f{G) corrections is 
negligible and all perturbation quantities except e follow 
their standard GR evolution). Note that this could be 
used as the initial condition for e when we evolve it nu- 
merically, and the initial condition for e could be obtained 
simply by taking its time derivative. 



The presence of the term 



27 1- 



16 e^Fo 



^ in Eq. (gl]) is no- 



This is not the whole story. We could use Eq. ^ 

to substitute the term -3^|Z - f^^ry in Eq. gB]) to 
re-express the evolution equation of e as 



e + 
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S' 



(47) 



where 
S' 



3k , f ( F F \ P 

^^"^'-'^''^-'[¥-39) 

9 ^ / 2F f\ k _^ 



The term (^1 + |l) in Eq. (gZl) makes the situation 

more complicated. During the whole matter-dominated 
era and part of the deceleration-to-acceleration transition 

period we have (l + |l) < 0. Thus, for the evolution of e 

to be stable on small scales [k"^ /o^H"^ ^ 1), we must also 
require that this term is subdominant compared with the 
third term in the squared brackets; that is, \FcH^ \ must 
be close enough to zero not only at early times, hut also 
at low redshifts {e.g., z < 0(10)). For example, in the 
ACDM hmit \FgH^\ = 0, so that the small-scale insta- 
bilities will never appear. In general, since the deviation 
from ACDM in the /(G) model is roughly characterized 
by the deviation of A from zero, A^ \ should be satisfied 
in order that the model evades cosmological constraints 



from linear spectra. 

To illustrate the effects discussed above, we have shown 
the evolution of some linear perturbation variables of the 
present model in Figs. [3] and ID Plotted in Fig. [3] are the 
evolutions of the (large-scale) Weyl potential 0fe. This 
is the coefhcient of the harmonic expansion of £ab as 
£ab = ^^k^^'t'kQ^b/a^ and is related to the Newto- 
nian potential, ^, by 4" = (f) — kRo?' / 2k^ for any specified 
/c-mode, where 11 is the anisotropic stress. From Eqs. ([TTl 
[33l- [36|) . it is obvious that depends on e, and so from 
the analysis above it is easy to understand why on smaller 
scales 4>k changes so dramatically. This situation is quite 
similar to that in the Palatini f{R) gravity model (l8| 
where f{R) = R + a{-R)^ with /3 > 0. Since the time 
evolution of (pk determines the CMB power through the 
integrated Sachs- Wolfe (ISW) effect as 



7-ISW 



where ji (fcr) are the spherical Bessel functions, and tq the 
conformal time at present, the extremely rapid variations 
in (pk might greatly enhance the angular power spectrum 
of temperature anisotropics, as in the the /3 > case of 
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FIG. 3: (Color Online) Evolution of Weyl potential ver- 
sus the cosmic scale factor a at different scales, k — 
0.001,0.002,0.005,0.01 Mpc"^ respectively from bottom to 
top . The values of A are indicated beside the curves. A — 
corresponds to the ACDM model. 




a 



FIG. 4: (Color Online) Evolution versus cosmic scale 
factor a of the cold dark matter density contrast 
AcDM at different scales, from bottom to top k = 
0.005,0.01,0.02,0.05,0.1 Mpc"^ respectively. The values of 
A are indicated beside the curves. ^4 = corresponds to the 
ACDM model. Note that the rapid growths of small-scale 
perturbations may make these scales leave the linear regime 
much earlier than in the standard ACDM model. 



Ref. [l^. This will provide the first stringent constraint 
on the present model (or equally on A). 

In Fig. m we have displayed the time evolution of the 
cold dark matter density contrast, Acdm, on different 
scales. As expected, on large scales the fc^ term in 
Eq. (|T7)l never becomes important and FgH^ is small 
enough for our choices of A, so the deviation from ACDM 
{A = 0) is small. On small scales, however, the term 
is significant, and makes Acdm blow up quickly. This is 
also similar to the /3 > branch in Ref. [3, and pro- 
duces a matter power spectrum which is strongly scale 
dependent for large k (note the difference to metric f{R) 
models [ill). This scale dependence will be strongly dis- 
favored by current data on galaxy power spectra (such 
as those from Sloan Digital Sky Survey) and gives a sec- 
ond stringent constraint on the /(G) cosmological mod- 
els. Considering these, although evaluating the numerical 
constraint on A is beyond the scope of the present work, 
we can claim that the parameter space for a viable f{G) 
cosmology is highly limited. This, among others, once 
more reveals the difficulties appearing in explaining the 
cosmic acceleration with modified gravity models. 

One should notice that the instabilities of the mat- 
ter component found in these /(G), and in f{R) models 
within the Palatini approach, are in fact dependent on 
the nature of the dark matter and not only on the grav- 
ity sector [l7j . For instance, instabilities are also found in 
GR models where cold dark matter is coupled to a light 
scalar field \3f\. However, when dark matter is not cold, 
and so has a free-streaming length, instabilities might 
go away, as occurs in some interacting hot-dark-matter- 
dark-energy models [ssj . Hence, ruling out these gravity 
models just due to the instabilities in CDM might not be 
the last word. 



V. CONCLUSION 

To summarize: in this work we have consider the cos- 
mology arising from a new modified gravity model, the 
modified Gauss-Bonnet model, where a function, /(G), 
of the Gauss-Bonnet invariant is added to the Einstein- 
Hilbert gravitational action to account for the current 
cosmic acceleration, at both the background and first- 
order perturbation levels. 

For the background- level evolution, we find that /(G) 
models cannot describe arbitrarily parameterized cos- 
mic histories leading to the current observations, because 
these histories are generally characterized by a transition 
from G > to G < which might impose too strong a 
requirement on the form of /. 

For the linear-level evolution, we present the first-order 
perturbation equations for the /(G) gravity model in the 
CGI formalism with a derivation which can also be ap- 
plied to general f {R,R'''' Rat, R"'""^ Rated) models. The 
special combination of R, R°'^Rab and R'^^'^'^ Rated terms 
in the invariant G ensures that the set of perturbation 
equations resemble that in the f{R) gravity models and is 
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much simpler than those in f{R°'''Rab) or f {R"'""^ Rabcd) 
theories. We analyze the perturbation equations and 
find that only the d^f/dG^ < subclass of the mod- 
els could have stable perturbation growth. Furthermore, 
even within this subclass, there will be a period during 
which the small-scale perturbation growth is no longer 
stable unless \Fa\H^ is close enough to zero. This un- 
stable growing period has two important consequences. 
Firstly, it makes the gravitational potential (pk change 
very rapidly, which may greatly enhance the ISW ef- 
fect and alter the CMB power. Secondly, the small- 
scale dark matter density perturbations grow much more 
quickly than in the ACDM paradigm, which might lead 
to a strongly scale-dependent matter power spectrum. 
In both cases, the model therefore faces stringent con- 
straints from the cosmological data sets on linear spectra. 
We thus conclude that the parameter space for a viable 
/(G) cosmological model is highly constrained. 

There are several points which arise from this work 
but are beyond the scope of the present investigation. 
For example, it might be interesting to make an analysis 
of the evolution dynamics for general /(G) models and 
to include the constraints arising from solar-system tests 



of gravity on the model. 
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